Time behaviour near to spectral singularities 
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Spectral singularities such as exceptional points invoke specific physical effects. The present 
paper focuses upon the time dependent solutions of the Schrodinger equation. In a simple model it 
is demonstrated that - depending on initial conditions - within close proximity of exceptional points 
the time behaviour of the wave function displays characteristic features such as very fast decay or 
the opposite, i.e. very long life time. At the exceptional point the wave function typically has a 
linear term in time besides the usual exponential behaviour. 
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INTRODUCTION 



Recently the effects of spectral singularities upon scat- 
tering have been discussed in general terms [l| using a 
particular model for demonstration. It is the exceptional 
points (EP) that generically occur in Hamiltonians of 
open systems for specific (complex) values of some pa- 
rameters 0, The EPs are associated with the co- 
alescence of two (or more) eigenvalues under variation 
of appropriate parameters; in contrast to a degeneracy 
the corresponding eigenstates also coalesce and have zero 
norm. For that reason they give rise to a double pole 
in the scattering function |5[ and thus to characteristic 
physical effects. Here we mention atomic and molecular 
physics optics jjj, nuclear physics [k| and in dif- 

ferent theoretical context 'PT-symmetric models [ll|, to 
name just a few. Depending on the particular situation 
EPs can signal a phase transition 



13(. The topo- 



logical structure being a square root branch point in the 
complex plane has been shown experimentally to be a 
physical reality HUB]. 



The present paper resumes the model dealt with in 
however, from a very different viewpoint. While the 
scattering deals with a stationary physical situation, we 
here investigate the effect of EPs upon the time behaviour 
of solutions of the time dependent Schrodinger equation. 
Such effects have been seen recently in [13] with optical 
micro-spirals, in 18[ for Rabi oscillations and earlier in 
19] for a singular di-electric tensor. It is expected that 
the rapidly improving experimental techniques of lasers 
applied for instance to atomic and molecular physics will 
focus upon time behaviour to an increasing extent [2cj |. 



The following section presents the time evolution with 
the emphasis at or close to an EP. A generic physical ex- 
ample is given in section 3 showing the characteristic and 
dramatic changes of the time dependent wave function in 
the proximity of EPs. A summary concludes the paper. 



TIME EVOLUTION 



The Model 



It is well established that, in the close vicinity of an 
EP, a two-dimensional matrix model suffices to capture 
all essential features associated with the singularity. We 
thus begin with the model Hamiltonian 



H{\) = Ho + Hi (A) = Hq + XV 
uii \ . / ei 6 



uj 2 I U £2 



(1) 



where the parameters ojk and determine the non-inter- 
acting resonance energies Ek = oJk + Ae^, fc = 1,2. The 
time evolution of a two dimensional state vector \tp{t)) 
with the initial condition \tp(0)) = {Ci, C 2 } T is given by 



\m) 



J2 eM-iEkWt){4>k\i>(0)}\4>k) (2) 



fe=l,2 

with the eigenvalues 

£?i, 2 (A) = -(wi + Wa + Afei + eijTD 

D = ^CC{\- EP\){\~ EP2)) 
CC = 46 2 + (e 1 -e 2 ) 2 

expressed in terms of the EPs 

i(ujx - L0 2 ) 



EP1 



EP2 = 



-28-i{e x -e 2 ) 

i(u)i - uj 2 ) 

e 2 ) 



+26 - i(ei 
and the normalised eigenvectors 
D 

UJ 2 



(3) 

(4) 
(5) 



(G) 
(7) 



2XS, 

A(ei - e 2 
1 



V>i - uj 2 + A(ei - e 2 )) 2 + (2A<J + D f 



(8) 



with \<j) 2 ) obtained from |</>i) by the replacement 8 — > —6. 
Note that the denominator of \<fik} vanishes when A — > 



2 



EP1 or A -> £P2, the leading order being (A - EP) 1 / 4 
1- 

Actually, as in [lj, we have chosen for the eigenvec- 
tors a basis rotated by the angle 7r/4 rather than the 
basis given by ([T]). This specific observational basis is 
essential to ensure that even for the non-interacting case 
(6 = 0) the two rotated channels feature equally and are 
observed simultaneously, while - for 6 = - the states in 
the basis given by ((T|) live in orthogonal spaces. The dif- 
ference between the two rotated channels appears when 
d is switched on as is discussed in the following sections. 

Inserting the explicit expressions into ([2]) we obtain for 
the components denoted by zi t2 (t) of the two component 
solution 



the result 

z 2 (t) = C 2 



( exp{—iEit) + cxp(— iE 2 t) 



Ci 



V 2 

exp(— iE\t) — exp(— iE 2 t) 

— 

D 

exp(—iEit) — exp(—iE 2 t) 
2D 



x (cji - lj 2 + A(ei - e 2 )). 



(9) 



The first component z\ (t) is obtained from z 2 (t) by the 
replacements C\ «-+ C 2 and S — ► —S. 

Before discussing the detailed behaviour in an actual 
physical situation we first turn to the analytic solution 
at an EP. 



Time behaviour at an exceptional point 

When A approaches an EP the denominator D in 
^ vanishes. So do the corresponding numerators since 
Ei — > E 2 when A — > EP. A finite limit is obtained and 
reads 

(ex — e 2 — 2iS + iS(uJi — 0J 2 )t)C 2 — 8{ui — uj 2 )tC± 



x exp(— it 



€\ — e 2 — 2iS 
i(eiuj 2 - t 2 jjJ\) + S(uji+ lu 2 ) . 
i{e 1 -e 2 )+2S ' 



(10) 



and 



(ei — e 2 — 2iS — i5(u)i — cu 2 )t)Ci — S(u)i — cu 2 )tC 2 



x exp(— it 



ei — e 2 — 2iS 
i(eiu 2 - t 2 io\) + (5(o;i + uj 2 ) , 
i(ei - e 2 ) + 25 ' 



The essential point is the linear time dependence oc- 
curring in the coefficient of the exponentials. The 'sin- 
gularity', usually associated with an EP, invokes in the 
time frame the additional linear time dependence. This 
is of course no surprise; it is a well known result from the 
theory of ordinary differential equations. The evolution 
equation idt\tp(t)) — H\ip(t) gives rise to linear terms in 
time if H cannot be diagonalised but has a non-diagonal 
Jordan normal form; this is precisely the case at an EP. 
It is here where the associate vector defined by 



(H(X = EP)-E EP )\ 



:) = \4>EP) 



(11) 



plays its role (see for instance [2l| ) . For more details and 
generalisations (higher order terms in t) we refer to the 
Appendix. 

We note that the linear time dependence at an EP has 
been noticed in for a non-diagonalisable di-electric 
tensor, and recently in fl7l ] for asymmetric scattering in 
optical micro-spirals, and also indirectly in lSJ for Rabi 
oscillations in a microwave cavity. 

We mention that these results cannot be obtained eas- 
ily from a simple Fourier transform of the corresponding 
Green's function G(E) = (E — H)^ 1 as given in [lj. The 
reason is due to the fact that (i) the residues blow up 
owing to the vanishing norm of the eigenfunctions at the 
EP (actually G{E) has a pole of second order Q), and 
(ii) the associate vector does not naturally occur. Of 
course, the mathematical connection does exist, yet the 
results in the time domain are much easier obtained by 
the approach used in the present paper. 



A PARTICULAR PHYSICAL CASE 

For illustration we use the same parameters as in [lj], 
that is ujx = 1.55-0.007i, u> 2 = l.l-0.007i, e 1 = -0.4- 
0.0006i, e 2 = 0.4 + 0.0005i and 5 = Q.0115*. For A we 
restrict ourselves to the interval [0.53,0.59]; the critical 
value is A c = 0.563, it is there where a resonance pole 
comes closest to the real axis. The close proximity of two 
exceptional points makes the movement of the resonance 
pole under variation of A rather swift and dramatic. This 
is illustrated in Figure 1 where trajectories are drawn in 
the complex energy plane. 

The specific values chosen have no particular signifi- 
cance, they serve to illustrate the principle, that is the 
effect of a near EP upon the behaviour of resonance poles 
and thus the time behaviour of the state vectors; any 
other set that invokes repulsion and coalescence of res- 
onances would serve the same purpose as long as the 
imaginary parts of the poles are near to the real axis. 
While wc did not specify units for the present illustra- 
tion, choosing for instance meV for the energies would 
imply picoseconds as time units in Figures 2 and 3. 

Similar in spirit to the procedure in 17[ we choose 
for the initial condition of ^(f) a specific polarisation 
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FIG. 1: Energy eigenvalue trajectories when A sweeps from 
0.53 to 0.59. Note that for our choice of parameters the tra- 
jectories run in opposite directions, the respective starting 
points are indicated by a dot. 



and consider the two cases ip(0) = {1,0} T and ij)(0) = 
{0, 1} T . In Figure 2 the two components are plotted for 
the latter case at a value of A where the two resonances 
are well separated (the dots in Figure 1). Three main 
features become obvious: (i) the initial conditions re- 
quire that the first and second component begin at zero 
and unity, respectively; (ii) the beat is determined by the 
frequency difference of the two near resonances while the 
coupling invokes an amplitude transfer between the two 
components; (hi) the damping - not to be confused with 
the beat - is determined by the width (the imaginary 
part) of the two resonances, exp(— tT) gives the envelope 
of the damped oscillation. Note that for our choice of 
parameters we find T = 0.007 for the width while the 
beat frequency AE = U[Ei] - R[E 2 ] = 0.025 is larger 
than the width (a few beats can be accommodated un- 
der the envelope) . Swapping the initial conditions would 
essentially swap the two components. 

A dramatic change occurs when is plotted at the 
critical point A c = 0.563 being the point where the tra- 
jectories peak in Figure 1. Moreover, the change affects 
the two different initial conditions in a very different and 
characteristic way. Three important changes are noticed 
in Figure 3: (i) the beat has disappeared; (ii) for the ini- 
tial condition i/j(0) = {0, 1} T (top row) the leading com- 
ponent is very weakly damped while the other component 
remains rather small, yet it is also weakly damped; (iii) 
for the initial condition ^(0) = {1, 0} T (bottom row) the 
leading component z\{t) is strongly damped while the 
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FIG. 2: Time dependence of the two components of the wave 
function at A = 0.53. The corresponding imaginary parts look 
the same. 



other component Z2(t) is essentially like that for the top 
row (here z\(t)). Again we emphasise that the respective 
envelopes of z 2 (t) (top) and Z\ (t) (bottom) are given by 
the damping related to the widths of the top and bot- 
tom peak in Figure 1, the values are r top = 0.0005 and 
r bot = 0.013. 

These findings are explained using Figure 1 and (I9I10[) 
and (jlip . At the critical point there is one narrow reso- 
nance (top peak in Figure 1) and one broader resonance 
(bottom peak in Figure 1). As the frequencies (the real 
parts of the energy) are the same, there is no beat. Now 
we first turn to z 2 (t) of the top row and Zi(t) of the 
bottom row in Figure 3. As seen from @ the expres- 
sion associated with C 2 = 1 has two terms: the sum of 
the exponentials with the two eigenenergies (and the fac- 
tor 1/2) and their difference (with the factor S/D); the 
two terms are added. Similarly, the corresponding two 
terms for z\ (t) being associated with Ci = 1, C 2 = are 
subtracted. The conspiring behaviour due to the sum 
in the one case and the difference in the other becomes 
evident in ([TOj) (with C 2 = 1, Ci = 0) and JTT]) (with 
G\ = 1, C 2 = 0), respectively. The energies in the ex- 
ponential are of course different in (|9]), there is the small 
and large damping (width). When adding the two terms 
as for z 2 (t), top row of Figure 3, the small damping pre- 
vails while subtracting as for Z\{t), bottom row, the large 
damping is dominant. At the top row Zi(t) is immaterial 
in comparison with z 2 {t), while at the bottom row z 2 (t) 
becomes stronger than Z\(f). We return to this interest- 
ing aspect in the last section. 

These rather dramatic differences for the time beha- 
viour depending on the initial conditions are expected to 
be observable as a physical effect; in fact, the polarisation 
in [17| is just one case in point. The findings discussed 
in the present paper complement beautifully the results 
presented in [l| . This is discussed in greater detail in the 
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FIG. 3: Time dependence of the two components of the wave 
function at the critical points A = 0.563. The top row is with 
initial conditions as in Figure 2 while in the second row the 
initial components are interchanged. Note that the time axis 
extends only half as far as the one in Figure 2. 
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following section. 

When going beyond the critical point with A, say at 
A = 0.59, Figure 2 is essentially repeated. The graphs for 
the imaginary parts of the complex components Zi t2 {t) 
look essentially as the ones given. 



SUMMARY AND CONCLUSION 

In [l[ the effect of EPs has been discussed within the 
same model and parameters as in the present paper, but 
the emphasis was on resonance scattering, that is a sta- 
tionary physical situation. Now we focus upon the re- 
lated time dependent problem, that is the solution of the 
time dependent Schrodinger equation, which is formally 
given by 

|^)>=exp(-ifft)|^(0)). 

For the purpose of the present paper the Hamiltonian is 
modelled by a 2 x 2-matrix describing an open system. It 
depends on a parameter A and of interest is the behaviour 
of \ip(t)) when A is at an EP or sweeping over a range 
nearby an EP. Similar to the findings in |l| the variation 
of A causes dramatic changes due to the proximity of 
EPs. 

When the resonances are well separated, the time be- 
haviour is damped owing to the widths of the eigenstates, 
but there is also a beat owing to the difference of the two 
frequencies. This changes dramatically when a pair of 
EPs is approached: the beat disappears, and depending 
on the initial condition the leading component is very 
weakly or very strongly damped. This is reminiscent of 
the very sharp resonance peak in the one channel and 
a broad peak in the other channel in two channel scat- 
tering [lj]. In addition, while the smaller component re- 
mains insignificant for the weakly damped case, it be- 
comes the stronger component for the strongly damped 
case. In fact, while there the leading component tends 
quickly to zero due to the larger width, the associated 
smaller component becomes dominant and is actually 
weakly damped. This corresponds exactly to the peak 
on the broad peak as detected in [l[ . 

As a general result we noticed that the singular beha- 
viour usually associated with EPs leads to a much milder 
pattern in the time frame. The 'self-orthogonality' of the 
eigenstates appears in a less dramatic way: depending 
on the order of the EP a polynomial in the time variable 
occurs besides the exponentials; in the simplest case it is 
just a linear term 17- 19j. 

As an increased activity of experiments using time be- 
haviour in atomic and molecular physics can be expected, 
the findings of the present paper may be of special rele- 
vance. 



Appendix 

This is a rehash of known facts from linear ordinary 
differential equations. If the time independent square 
matrix operator O is diagonalisable, the evolution equa- 
tion 



o\x(t)) 

with the initial condition |x(0)) is solved by 

>llx(o)M) 



lx(*)> 



k 



exp(E k t)- 



with Ek and \4> r k ) and (<p k \ forming the right and left 
hand eigensystem of the operator O. At an EP, when two 
eigenvalues coalesce, the matrix cannot be diagonalised 
but has the Jordan decomposition 



o = sjs- 



with 



(E EP 1 
E EP 



En 



V 



E N ) 



where we have listed the two EPs first. Now the solution 
of 



reads for the initial vector {Ci, . . . , Cjv} t 



|£(t))= { {Ci+tC 2 )exv{E E pt),C 2 exp(E EP t), 
C 3 exp{E 3 t), ...,C N exp{E N t)} T 

which can be transformed back into the basis of the orig- 
inal O using the similarity transformation S. A so-called 
associate vector related to the eigenvector \4>ep) by 



>EP 



features here in the second column of S. For the special 
case of Section 2 it can be chosen as in the second column 
of S which reads 

/■ 2<5+z( ei -e 2 ) \ 

S= 

V ) 

This explicit form explains the terms linear in t occurring 
in (HI and (fTTjl . 

It is clear how to generalise this when more than two 
levels coalesce. For k coalescing levels the first compo- 
nent of is given by 

(Ci + 1 C 2 + t 2 /2l C 3 + ... + t^/ik - 1)! C k ) exp{E EP t) 



■5 



and accordingly for the following components. 
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